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Abstract. We fix a prime number p and K a number field, we denote by 
M the maximal abelian p-extension of K unramified outside p. The aim of 
this paper is to study the Z p -module X = Gal(M/i^) and to give a method 
to effectively compute its structure as a Z p -module. Then we give numerical 
results, for real quadratic fields, cubic fields and quintic fields, together with 
interpretations via Cohen-Lenstra's heuristics. 

1 Introduction 

We fix a prime number p and a number field K . We denote by M the maximal 
abelian p-cxtcnsion of K unramified outside p. The aim of this paper is to 
study the Z p -module X = Gdl(M / K) and give an algorithm to compute its 
Zp-structure. This module is described by the following exact sequence from 
class field theory ( [Grail p. 294]) 

Uk -n„|pEtf »Gb1(H/K) (1) 

where Uk is the pro-p-completion of the group of units Uk, U„ is the group 
of principal units at the place v above p of K, and Ti is the maximal p-sub- 
extension of the Hilbert class field of K. Leopold't conjecture for K and p is 
equivalent to injectivity of Uk Yl v \ P Uy- Therefore, from this exact sequence, 
we deduce that the Z p -rank r of X is greater or equal to r-i + 1 and is equal 
f'2 + 1 if and only if K and p satisfy Leopoldt's conjecture. Hence X is the direct 
product of a free part isomorphic to Z^ and of a torsion part, that we denote 
by 7p. Our algorithm checks if K satify Leopoldt's conjecture at p and then 
compute the torsion T v . 

We propose a method which is based on the fact that the Z p -module X is 
the projective limit of the p-parts of the ray class groups modulo p n , A p ™(K). 
We then study the stabilization of these groups with respect to n and the be- 
haviour of invariants of A p n(K), as n is increasing. This approach leads us to 
our algorithm. 

Before addressing the technical part of this article, we recall the definition 
and some basic properties of the ray class groups modulo p n . Then, we use our 
algorithm to compute some cases and propose an heuristic explanation of the 
statistical data, using the Cohen-Lenstra philosophy f |C-L| ). 
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2 Background from class field theory ( fGral] , [Ser] ) 



In this section, we recall the basic notions from class field theory that we will 
need later. We fix v a place of K above p and ir v a local uniformiser of K V) the 
completion of K at v. 

Definition 2.1. 

1. The conductor of an abelian extension of local fields L v /K v is the min- 
imum of integers c such that U° C N Lv / Kv (L*) (we recall that U£ 

1 + (ny) and we use the convention U® — U v ). 

2. The conductor of an abelian extension L/K of a global field is the ideal 
m = Y[ v Pv v > where v runs through all finite places of K and where c v is 
the conductor of the local extension L„ / K„. f\Grall p. 126-127] Theorem 
and Definition 4-1 + Lemma 4.2.1). 

We start with 2 lemmas. 

Lemma 2.2. ('Ser! p. 219). Let K v be the completion of K at the valuation v 
normalized by v(p) = 1 and v(ir v ) — where e v is the ramification index of the 
extension K v /Q p . If m> -^zj, then the application x —¥ x p is an isomorphism 
from U™ to U r v n+e v . 

Lemma 2.3. Let K v C L v C M v be a tower of extensions of Q p , such that 
the extension M v /K v is abelian and the extension M v /L v is of degree p. We 
denote respectively by cm,v and cl, v the conductors of the extensions M v /K v 
and L v /K v . If ci,,v > -^zj, then we have 

cm,v < cl.v + e„. 

Proof. By definition cl, v is the smallest integer n such that U™ C Nl v /k v (L*)- 
Local class field theory gives the following diagram 

1 N Mv/Kv {M*) »- K* - G&\{M V /K V ) »- 1 

1 >- N Lv/Kv (L*) K* > G&\(L V /K V ) 1 

Applying snake lemma we get the exact sequence 

1 N Mv/Kv (M* ) ^ N Lv/Kv (L* ) ^ G&\(M V /L V ) = Z/pZ > 1. 

Consequently N Mu / Kv (M*) is a subgroup of index p of N Lv / Kv (L*). Let 
n € N, n > cl,v + e v and x € {/". We have to show that x € N Mv / Kv (M*). By 

— _ 
lemma 12721 xp is a well defined element of G f7™ e ". Yet n — e v > cl, v therefore 

xp € N Lv / Kv {L*). Now, as N Mv / Kv (M*) is of index p in N Lv / Kv (L%), we 
deduce that x £ N Mv/Kv (M*). We have therefore Uff C N Mv/Kv (M*) for all 
integers n such that n > cl + e v . By definition of the conductor, this proves 
cm.v < cl, v + e v . 

□ 
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Definition 2.4. Let n be a positive integer. We denote by 

• H the maximal abelian unramified extension of K ; 

• H p n the compositum of all abelian extensions of K whose conductors divide 

P n ; 

• "Hp" the compositum of all abelian p-extensions of K whose conductors 
divide p n ; 

• M the maximal extension of K which is abelian and unramified outside p. 

So the Galois groups G&l(T-L/K) and Gai(H p n/K) are respectively isomorphic 
to the p-parts of G&\{H / K) and G&\{H p n/K). 

Proposition 2.5. (\ Grail p. J^l] corollary 5.1.1) We have the following exact 
sequences, 

1 _ K < U v IL b U™« 1 K Gal(JW*) — - 1 

1 \[ v Uv Tk G&\(H/K) 1, 

where Xk is the group of ideles of K . 

We denote the Galois group Gal("Hp» / K) by A p ^ (K). It is the p-part of the 
Galois group Gal(_ff p n / K) which, in turn, is isomorphic to the ray class group 
modulo p n of K. By definition, we have a natural inclusion H p ™ C TL p ™+i , the 
union (^J'Hp™ is equal to M and the projective limit hmApn(K) is canonically 

n 

isomorphic to X. 

Proposition 2.6. For any integer n > 0, the Galois groups of the extensions 
M and H p n of K are related by the following exact sequence 

1 UP U V \ P Gal(M/X) Gal(H p n/K) 1, 

where U£ ^ = {u 6 Uk such that Vv\p, u 6 U™ £v } and 

uf ] — - ELip X A p n (K) 1, 

where U% "* is the pro-p- completion of \ i.e fim ^ /p m . If Moreover K 

m 

( p n ) 

and p satisfy Leopoldt' conjecture, then U K — > Ylv\p ^" e " * s injective. 
Proof. To obtain the second exact sequence, we apply pro-p-completion pro- 
cess to the first. Note that injectivity of U K —> Yl v i p U" Sv is equivalent to 
Leopoldt's conjecture. Now we prove exactness of the first sequence. 
From the definition of the extensions M and H p n , we deduce the diagram 

1 K* U v n Hp 1 Ik Gal(M/X) 1 



K x EU U v EU U: e * Ik Gal(^/X) ■ 
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It follows immediately from the snake lemma that 



ker(Gal(M/X) -+ Gal(H p »/K)) = (K x J] U v [J U^)/(K X J] ^ II 
Now, we define the application 

o: (K* n u v n ^" e ") -»• (II wrvuTl 

v\p v\p v\p 

by setting for fc(u„)„ G K x ]\ V ^U V ]\ V ^U^ , 6(k{u v ) v ) = {u v ) v \ p , where 
(u v ) v \ p is the class of (u„)„| p in (n„| p U™»)/U$ } . 

We first check that the application is well defined, i.e. that if k(u v ) v = 
k'(u' v ) v in K x rL +P ^IL b £C% then 9{k{u v ) v ) = 9(k'(u' v ) v ). By definition, 
for all v, k(u v ) v — k'(u' v ) v if and only if i v (k)u v — i v (k')u' v , where i v is the 
embedding of K in K v . We deduce that for all v, i v (k'k^ 1 ) e U v and that for 
all v\p, € t/" e ". So we get /c'fc -1 e fT^ 5 ) and (u„)„| p = «)t,| p - 

It is clear that (K x Iltffp^f IIw|p 1) c ker(#) and that the application 9 is 
surjective. We will show that (K x Yl^ p U v Yl v \ p 1) = her (6). Let k(u v ) <G ker(#), 

there exists x € ' such that for all t>|p, tt„ = i v (x). We consider the 
element x(u' v ) v , where u' v = 1 if u|p and u^, = i„(a;) u„ if u j p. We have 
(u v ) v = x(u' v ) v fc(u„)„ = fcx«)„ and as kx(u' v ) v e {K x \[ vl[p U v ]J vlp 1), we 
have kcr(9) C (i^ x ]X tp ^ lL| P 1) and finally 

(k x n u v n ^r-)/(^ x n ^ n - n ^i^p. 

v\p v\p v\p v\p v\p 

The result follows. □ 



3 Explicit Computation of T p 

In this section, we present our method to check that K verify Leopoldt's con- 
jectue at p and then to compute T p . The main point is that, for n large enough, 
Apn{K) determines X. 

3.1 Stabilization of A p n(K) 

For simplicity we note Y n — ker(_4 pn +i (K) — > A p ™(K)). Let K be the com- 
positum of all the Z p -extensions of K. We denote by r the Z p -rank de X, so 
r > r 2 + 1. 

Proposition 3.1. There exists n such that K n V. p ^a /K r\H p is ramified at 
all places above of p and for all n ^ n , Y n surjects on (Z/pZ) r . 

Before proving the proposition we need a lemma. 
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Lemma 3.2. If the extension K (~l % p n / K n H p is ramified at a place v above 
p, then c n , v > where c„ tV is the conductor of the local extension [K f~l 

7~Lp™ )u; j Ky ■ 

Proof of the Lemma \3.2\ As M contains the cyclotomic Zp-extension, there ex- 
ists a no such that K n H p n /K fl % p is ramified at all places u above of p. As 
KCiHp'io /Kr\W p is ramified at w then, for n > n , K n"H p n / K CiTip is ramified at 
v, so that there exists to such that n > m > 2 and that KnW pm -i jKCvhip is un- 
ramified at u and such that KfXHpm /KfXHp is ramified at v. Then, the local con- 
ductor c„ hV is greater than (to— l)e v , yet to > 2 so c„ ltV > (m—l)e v ^ e v ^ 

As the conductor of the local extension K n / if divides the conductor of 
K n Hpn/K, we have c„,„ > cm, w > □ 

Proof of the Proposition \3.1l We consider the following diagram. 

K n H p « (K n H p n )H P H p n (2) 




K 

We have Gal{K/K) = Z£. It is clear that Y n -» Ga\(K n H p n + i/K n 
% p »). Yet Gal(if/if n ?V) is a Z p -sub-module of Gal(if/iT) = Z p of finite 
index, therefore it is isomorphic to Z p . Hence there exists r extensions, say 
Mi, M 2 , • • • , M r of if n H p », contained in if such that Gal(M;/if H % p ») ~ 
Z/pZ and Gal(M 1 ---M r /K n H p ») ~ (Z/pZ) r . Yet the conductor of the 
extension if n % p n/K divides p n = n u | P P™ e "- Moreover the hypothesis on 
if n "H P " /K n H p ensures that we can use Lemma 12.31 and consequently the 
conductor of the extension Mi/K divides Ylvlp P™ e,,+e " = P™ +1 j i- e - Mi C H p *.+i 
for all i G {1, • • • , r}. Hence the map is surjective. □ 

We deduce immediately the corollary. 

Corollary 3.3. We assume that for a naturel number n the extension if fl 
Hpn /if n H p is ramified at all places above of p, and that the cardinal of Y n is 
exactly p T2+1 . Then Y n ~ (Z/pZ) r2+1 and if verify the Leopoldt's conjecture at 
P- 

From now on, as we can numerically check that if sastify the Leopoldt's 
conjecture at p, we assume it, and we use it to compute T p - Note that if the 
Leopodt's conjecture is false, then r > r 2 + 1 and our algorithm never stops. 

Corollary 3.4. We assume that, for some integer n such that the extension 
if H H p n /K fl "Hp is ramified at all places above of p, the cardinal of Y n is 
exactly p r2+1 . Then, Y n ~ Gsd(K n H p »+i /if fl H p n) . 
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It remains to check that if Y no ~ (Z/pZ)'' 2+1 for some no, then Y n ~ 
(Z/pZ) r2+1 for all integer n > n . For this purpose, we consider the exact 
sequence defining the p-part of the ray class group. 

1 0, b tC 6 " ^ A p n (K) 1, 

and we note Q n = ]J v]p U™« /U { £ ) . We have Q n = G&1(M/H p n) and conse- 
quently Q n /Qn+i = Y n ~ Gal(H p n,+i/Hpn). 

Proposition 3.5. For n > 2, raising to the p th power induces, via the Artin 
map, a surjection from Y n to Y n+ \ . 

Proof. We recall that Q n = H v{p U™ e - /U { f ] = ker(X -> A p n(K)). We have 
that n > j^j. Raising to the p th power realizes an isomorphism of Y\v\ P U" ev 
onto Y\ v i p J/™ 6 " 4 " 6 ". This isomorphism induces a surjection from Q n onto Q n +i- 
We consider finally the diagram 



Q n+1 Q n Q n /Q 



(.)» 



(•)* 



n+l 



(.■)* 



1 >" Qn+2 >" Qn+1 >" Qn+l/Qn+2 >" 1 

We deduce from the snake lemma that the vertical arrow on the right side is a 
surjection from Q n /Q n +i onto Q n +i/ Qn+2, i-e. from Y"„ onto Y n+ i. □ 

Corollary 3.6. FFe denote q n = #(Y n ). For all n ^ 2, q n ^ (fo+i- Therefore 
the sequence (q n ) n >i is ultimately constant. 

We recall that Y n is kcr(.4 pn +i (K) -> ^n^)). 

Theorem 3.7. There exists an integer uq such that Y na ~ (Z/pZ) r2+1 . More- 
over for all integers n > n , the modules Q n — Gal(M/U p n) are Z p -free of rank 
r 2 + 1 and 

r„~ (z/ P zp +1 . 

Proof. The Z p -module X is isomorphic to the direct product of its torsion part 
and of Z p 2+1 . An isomorphism being chosen, we can identify Z p 2+1 with a 
subgroup of X and therefore define, via Galois theory, an extension M' of K 
such that G&\{M'/K) ~ 7^ and ifM' = M. 

This extension being unramified outside p, there exists an integer m such 
that M' C and consequently % p ^iK = M. Moreover for all integer 

n > m, Gal(M/U p n) is a sub-module of finite index of Gal(M/M') = Z£ 2+1 , 
consequently Q n = Gal{M/U p n) ~ Z£ 2+1 . The Z p -module Q n is therefore free 
of rank r 2 + 1. 

About the other kernel Y n we saw that there exists an integer n 2 such that Y n 
maps surjectively onto (Z/pZ) r2+1 for all integer n > n 2 (we can choose n 2 to 
be the minimum of all integers n such that for all p-places v the conductors 
of (K n H p n) w /K v are greater than or equal to -^). Finally we note that 
raising to the p th power realizes an isomorphism between U™ e " and f7™ e " +e ", 
hence the quotient Q n /Q n +i, which is isomorphic to Y n , is killed by p. Define 
no = Max(ni,n 2 ) and fix an integer n > no. The kernel Y n is therefore a 
quotient of Z p 2+1 , which maps surjectively onto (Z/pZ) r2+1 and is killed by p. 
Hence we get Y n ~ (Z/pZ) r2+1 . □ 
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3.2 Computing the invariants of T v 

We start by recalling the definition of invariant factors of an abelian group G. 

Definition 3.8. Let G a finite abelian group, there exists a unique sequence 
Oi, • • ■ ,at such that for all i, Oj|aj+i for i S {1, • • • ,4—1} and G ~ ]X' = i ^/fli^. 
These ai are the invariant factors of the group G. 

In what follows we will note them J-I(G) — [a±, ■ ■ ■ , a t ]. If G is a p-group, 
these invariant factors are all powers of p. In practice, we are able to determine 
the invariant factors of A p ™ (K) . We will see in this section that the knowledge of 
invariant factors of A p n(K), for n large enough, combined with the stabilizing 
properties of A p ™(K), does determine explicitly the invariant factors of, and 
thus T p - We recall that for n large enough, A p n(K) is isomorphic to the direct 
product of Gal(K n H p »/K) and of Gal(U p n/K n H p n) = T p . So we will first 
explore the structure of G&\(K n H p n/K). 

Proposition 3.9. Let uq be such that Knl-lp^o /KDHp is ramified at all places 
above of p and 

Y no ~ (Z/pZ) Ta+1 . 
Then for all integer n > no, we have 

G&\{K/K n H p n+i) =pGal(K/K C\U pn ). 

Proof. By Theorem 13.71 on one hand, Q n is Z p -free of rank r-i + 1 and on 
the other hand Y n = Q n /Q n+ i ~ (Z/pZ) r2+1 . This gives Q n+1 = pQ n . As 
Kr\W p no /Kfl'Hp is ramified at all places above of p and Y no ~ (Z/pZ) r2+1 , we 
have Tp C A P 'io (K), so KH p ^n = M. Then, considering the following diagram, 



K M 




K n H p n % 



K 

we get the required isomorphism. □ 

Corollary 3.10. Let no be an integer such that K fl "H P "o / K n Ti. p is ramified 
at all places above of p and such that Y no ~ (Z/pZ) r2+1 . Then for all integers 
n > no, i/ie invariant factors of Gal(K (~) H p n+i / K) are obtained by multiplying 
by p each invariant factor of Gal(K D T-L p n /K). 

From the fact that X ~ Zp 2+1 x T p , the ray class group, Gai(Hpn/K), is 
isomorphic to the direct product of Ga\(K n H p n /K ) and of G&\^H P ^ / K flHpn.). 
The invariant factors of Gal("H p n /K) are then simply obtained by concatenating 
those of the two groups forming the direct product. We now state the result 
that explicitly determines 7~ p - 
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Theorem 3.11. Let n such that Y n = (Z/pZ) r2+1 and K n H p n/K n H p is 
ramified at all places above of p. We assume that 

TT{A p n{K)) = ,h,ai,--- ,o r2 +i] 

with (up(oi)) > (up(6t)) + 1; 

J*X(.4pn+i (if )) = [6i,--- ,6 t ,pai,--- ,po r2+ i]. 

Then, we have 

JT(T P ) = [bi,--- ,6*]. 

Proof. Indeed, as 

r„~ (z/ P zr +1 , 

we have „4 p ; (if) ~ for Zp (X) x Gal(if n % p i/K) for ie{ii,n + 1}. We saw that 
the invariant factors of Gal (If n"H p n+i /if) are exactly equals to p times those of 
Gal(if nUpn /if). Consequently, if a is an invariant factor of Gal(if HUpn+i /if), 
we have necessarily a — pai or a — pbi . 

But as Min(up(aj)) > Ma,x(v p (bi)) + 1, none of the invariants factors of Gal(if n 
H p n+i/K) is of the form phi. The invariant factors of Gal(X n H p n+i/K) 
are therefore exactly pa\,--- ,pa r2+ \. The result follows from the fact that 
A p n+i (if) is isomorphic to the direct product of T p and of Gal(if fl H p n+i/K). 

□ 



4 Explicit computation of bounds 

More generally if we note e = max^ipje,,} the ramification index of if/Q and 
s the p-adic valuation of e, then we start to check whether A p ™ (if) stabilizes 
from rank n = 2 + s. To show that n — 2 + s is the proper starting point we 
consider the diagram, 




The places above of p are totally ramified in K s+ \/K s therefore Kr\'H p s /K(l 
H p is ramified at all places above of p and we start the computation checking 
whether A p n(K) stabilizes from n = s + 2. We first prove that all places above 
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p are totally ramified in K s+ \/K s . 
Considering the following diagram, 



K s+1 




Q 

the ramification index of p in Q s+ i/Q is p s+1 while the one in K/Q is p s a 
with p \ a. Therefore the extension K s+ i/K is ramified and K s+ i/K s is totally 
ramified at all places above p. 

Corollary 4.1. Let e be the ramification index of p in K/Q and s be the p-adic 
valuation of e. Let n > 2 + s, we assume that 

FI{Ap-n(K)) = [bi, ■■ ■ ,b t ,ai, ■ ■ ■ ,a r2+ i], 

with Min(u p (aj)) > Max(w p (6i)) + 1, and moreover that 

Fl(A p n+i(K)) = [bi, ■ ■ ■ ,b t ,pai, ■ ■ ■ ,pa r2+ i}. 

Then, we have 

T2(T P ) = [&!,-•• ,6*]. 

All the computations are been done using the PARl/GPsystem [PARI-GPJ. 

Example 4.2. We consider the field K = Q(v— 129) and p = 3. We have: 
T1(A P 2(K)) = [3,3,9], T1{A P 3{K)) = [3,9,27] and F1{A p a{K)) = [3,27,81]. 
We deduce that l~ p = (Z/3Z). 

5 Numerical results 

In the section, we give an explanation of some numerical results that we have 
computed. 

5.1 Heuristic approach 

We first recall some results on Cohen-Lenstra Heuristics. The main reference 
on the subject is the seminal paper of Cohen-Lenstra |C-L) , see also [Del . If we 
assume that the extension K/Q is Galois with A = G&\(K/Q), then the module 
T v is a Z[A]-module. In this section, we assume that A is cyclic of cardinality 
I, for some prime number Then T p is 0;-module, where Oi is the ring of 
integers of Q((i). In general, we know that all Oj-module G can be written in a 
non-canonical way as ©* =1 Oz/cii, where the are ideals of Oi. Yet the Fitting 
ideal o = Yii=i a i depends only of the isomorphism class of G, considered as a 
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0/-module. This invariant, denoted by a(G), can be considered as a generaliza- 
tion of the order of G. We also have N(a(G)) = #G. 

To simplify the notation we set. 



Eg n = Eg N(a(G))<N' wnere the sum is over all isomorphism classes of 
Oi-module G; 



• J2a,N — J2a,N(a)<N' 

• Ea'.AT ~ Ea' ,N(a')<N and a'Ap=l> 

• Ep n = Ep n(p)<n and pes ' where S p designed the set of all p-places of 
Oi. 

We consider a function g, dehned on the set of the isomorphism classes of 
0/-modules (typically g is a characteristic function). We then put 



S N (g) = E 
Sn = E 



g(G) 



G,N #Auto, (G) ' 



1 



G,JV #Aut 0( (G)- 

Definition 5.1. TTie average of g, if it exists is, the limit when N — > oo of the 
quotient 

Sn(9) 
Sn 

This average is denoted by Mi_o(g). 

As in |C-Lj . we denote by w(a) = J2g a(G)=a #Aut (G) ' wnere a is an ideal 
ofO,. 



Proposition 5.2. ( IC-LJj p. 40 corollary 3.8) Let neN, then 

The notation p Q ||a meaning p"|a and that p Q+1 { a. Consequently the function 
w, defined on the set of ideals ofOi, is multiplicative. 

Notation. We denote by n p the characteristic function of the set of the 
isomorphism classes of groups whose p-part is non trivial. 

Proposition 5.3. f\C-Lf p. 41 example 5.10) We note pi, •• • ,p ff the p -places 
of Oi, the average ofH p exist and we have 

i=l fc>i V * / 

where the fi designed the degree of the residual extensions Oi/pi over¥ p . 
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Corollary 5.4. If the extension K is a Galois extension, all residual degrees 
are equals to f and in this case 

\k>i f 

Remark. The real M; :0 (n p ) is called the 0-average. This notion can be gen- 
eralized to u-average. The expression to compute the u-average is obtained by 
replacing k by k + u in the expression of the 0-average. 

5.2 Somes numerical results 
5.2.1 Case of the quadratic fields 

We consider all quadratic fields of the type Q(Vd) with d square-free and < 
d < 10 9 . Then, we compute the proportion of fields with non trivial Z p -torsion of 
X. We note this proportion f exp . The relative error \f exp — M2 i o(n p )|/M2 i o(n p ) 
is denoted by S. 



p 


M 2 , (n p ) 


fexp 


5 


2 


0,71118 


0,93650 


0,31683 


3 


0,43987 


0,50120 


0,13942 


5 


0,23967 


0,23854 


0,00470 


7 


0,16320 


0,16280 


0,00247 


11 


0,09916 


0,09893 


0,00243 


13 


0,08284 


0,08266 


0,00212 


17 


0,06228 


0,06214 


0,00233 


19 


0,05540 


0,05526 


0,00260 


23 


0,04537 


0,04527 


0,00207 


29 


0,03375 


0,03560 


0,00193 


31 


0,03330 


0,03323 


0,00219 


37 


0,02776 


0,02770 


0,00198 


41 


0,02499 


0,02493 


0,00207 


43 


0,02380 


0,02376 


0,00152 


47 


0,02173 


0,02168 


0,00207 
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We consider now the quadratic field of the type Q(V") with — 10 9 < d ^ 0. 
One uses the 1-average that denoted by Af2,i(n p ). 



p 


M 2 ,i(n p ) 


fexp 


5 


2 


0,42235 


0,93650 


1.12734 


3 


0,15981 


0,25718 


0,60926 


5 


0,04958 


0,04909 


0,00989 


7 


0,02374 


0,02365 


0,00374 


11 


0,00908 


0,00905 


0,00416 


13 


0,00641 


0,00638 


0,00360 


17 


0,00368 


0,00365 


0,00445 


19 


0,00292 


0,00291 


0,00589 


23 


0,00198 


0,00197 


0,00510 


29 


0,00123 


0,00122 


0,00916 


31 


0,00108 


0,00107 


0,00929 


37 


0,00075 


0,00074 


0,00813 


41 


0,00061 


0,00060 


0,00982 


43 


0,00055 


0,00055 


0,00998 


47 


0,00046 


0,00046 


0,01626 



We have also computed these proportions for other fields and we consider 
the distribution of torsion modules with respect to invariants factors but they 
will not be given here. 



5.2.2 Explanation of numerical results 

In this section we explain our numerical result. Looking at two tables in ^5.2. II 
we remark that the proportion f exp for real quadratic fields seems to be a 0- 
average and a 1-average for the imaginary quadratic. We remark also that the 
default S for p — 2, 3 increases with the number of fields computed. To explain 
these phenomena we recall a computation of Gras |Gra2j p. 94-97. Let k be a 
number field, we denote by K = k(^ p ) and w the idempotent associated with 
the action of Gal(K/k) on 

Theorem 5.5 (Corollaire 1 p. 96 |Gra2j ) . Let p be a prime, p ^ 2. // \i v <f_ k 

then the torsion of X is trivial if and only if any prime ideal of k dividing p is 
totally split in K/k and {CIk) 1 ^ is trivial where CIk is the p-part of the class 
group of K . 

In case of quadratic fields, if p > 3 then \x p (jLk and the ramification index of 
p in Q(Cp)/Q is p — 1 ; then all prime ideal of k dividing p ramifies in K, therefore 
they are not totally split, hence the torsion is trivial if and only if (CIk) 1 ^ is 
trivial. So in case k is a real quadratic field the computation of T p reduces to 
the computation of a class group of imaginary quadratic field and we use the 
0-average following Cohen-Lenstra Heuristics. In case of imaginary quadratic 
the remark |Gra2j p. 96-97 explains the 1-average. In case p = 3, if d = 6 mod 9 
then the ideal of k above p is totally split in K, so the torsion is non trivial. It 
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explains why the frequency obtained is greater. If we consider the other average 
M^(n 3 ) = M 2)0 (n 3 ) x I + §; then we obtain 



N 


M^(n 3 ) 


fexp 


5 


lO 6 


0,50989 


0,48094 


0,05678 


10 v 


0,50989 


0,49054 


0,03794 


10* 


0,50989 


0,49697 


0,02533 


10 9 


0,50809 


0,50120 


0,01704 



We now make the computation without the d = 6 mod 9. 



N 


M 2 , (n 3 ) 


fexp 


5 


10 B 


0,43987 


0,40679 


0,07521 


10' 


0,43987 


0,41776 


0,05027 


10 8 


0,43987 


0,42511 


0,03356 


10 9 


0,43987 


0,42995 


0,02257 



It remains to study the 9-rank in case where d = 6 mod 9, and to try and 
find density formulas for the 9-rank. Finally the discrepancy in the case p = 2 
is explained by genus theory. Indeed if the discriminant is divided by enough 
primes then the torsion is not trivial, this explains that the frequency tends to 
1. 
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